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Abstract 

We establish the exponential time decay rate of smooth solutions of small amplitude 
to the Vlasov-Poisson-Fokker-Planck equations to the Maxwellian both in the whole space 
and in the periodic box via the uniform-in-time energy estimates and also the macroscopic 
equations. 



1 Introduction 

In this article, we study the convergence to the equilibrium for the nonlinear Vlasov-Poisson- 
Fokker-PIanck (VPFP) equations in the whole space. The VPFP system is one of the fun- 
damental kinetic models in plasma physics to describe the dynamics of charged particles 
(electrons and ions) subject to the electrostatic force coming from their Coulomb interaction 
and to a Brownian force modeling their collisions (Fokker-Planck). The VPFP system reads 

ft + V-V.J + diM{E - Hf) = DAJ (1.1) 

where / = f{t,x,v) > is the distribution of particles at time t, position x and velocity v 
for {t,x,v) G R4. X X M^, /3 > is the friction coefficient, -D > is the thermal diffusion 
coefficient, and E = E{t,x) is the self-consistent electric force. The equation (II. ip is coupled 
with the Poisson equation 

E = -V^, -A$= / fdv-1, ^ — >Oas|x|^oo (1.2) 

where $ = $(t,x) is the internal potential. 

The global equilibrium of (jl.ip and ()1.2p is given by the Maxwellian: 

2ttD 3 /3|i'i^ 

/ = ^:=^(^;):=(__)-5e-— , E = 0. (1.3) 

For simplicity, we take (3 = 1, D = 1. Letting 

f = ^x + ^/Jlg , 
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the VPFP system is written in the fohowing perturbed form: 



gt + v-V^g + v^-V^c/} + {-g - Vv9)-^x(i> - Lg = 0, 



A(f> 



g^dv, 



(1.4) 



where L is the hnearized Fokker-Planck operator given by 
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The Fokker-Planck operator is a well-known hypoelliptic operator. Diffusion in v together 
with transport v ■ Vx has a regularizing effect not only in v but also in t and x. Note that this 
is nontrivial since the diffusion only acts on the velocity. This phenomenon can be obtained 
by applying Hormander's commutator (cf. [15j) to linear Fokker-Planck operator. For more 
details, we refer to [ij. In particular, in smoothing property has been shown for the 
linear Vlasov-Fokker-Planck as a hypoelliptic operator when the external potential satisfies 
a certain condition. 

On the other hand, the Vlasov-Fokker-Planck (VFP) operator is also known as hypoco- 
ercive operator, which concerns rate of convergence to equilibrium. (See [23j). Indeed, the 
trend to the equilibrium with the rate of t"^/*^ is investigated in [10] in the case of the linear 
VFP equation with the external potential which is strictly convex at infinity. Later, it is 
proved in |16j that the solutions for the linear VFP equation with the external potential of 
high-degree approach exponentially to the Maxwellian equilibrium and the rate was given 
explicitly. In f23], Villani shows almost exponential decay to the equilibrium for the weakly 
self-consistent VFP equation in the periodic box with the Coulomb interaction potential re- 
placed by a small and smooth potential. When the coupling of the Poisson equation comes 
into play, then the fully nonlinear VPFP system adds more difficulties to study in this context. 
In [4], the convergence of free energy solutions in L^(M^) to the equilibrium for the VPFP 
equations with the confinement by the external potential was studied, but the convergence 
rate was not considered. 

To our knowledge, the convergence to the Maxwellian has not been yet shown for the fully 
nonlinear VPFP equations in the whole space without any confinement. The main goal of 
this paper is to establish the exponential convergence rate of the solutions for the nonlinear 
VPFP equations to the Maxwellian in the whole space in the regime of small and smooth 
solutions. 

Before we state the main results, we briefly review the existence theories for the VPFP 
equations. Global existence of the solutions to the VPFP system in two and three dimensions 
have been studied by many authors and we will not attempt to exhaust references in this 
paper. Classical solutions were obtained by Victory and O'Dwyer |25j in two dimensions 
and by Rein and Weckler 120) for small data in three dimensions. Bouchut [2] established 
the existence and uniqueness of a global smooth solution in setting in three dimensions. 
Asymptotic behaviors and time decay of the solutions near vacuum regime have been consid- 
ered by Caprio [5], Carrillo, Soler and Vazquez [8], and Ono and Strauss [19]. We mention the 
work of Victory [23], Carrillo and Soler [7J where the global weak solutions were constructed 
and also Bouchut [3] where the smoothing effect was observed. Recently, stability of the front 
under a Vlasov-Fokker-Planck dynamics was studied in |12] . 

There have been a lot of recent progress made on the results towards the convergence 
rate to the equilibrium near Maxwellian regime for other fundamental kinetic models such as 
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the Boltzmann equations with various colhsion kernels, for instance see [9l [T8| I22|. We 
also mention that there are quite recent results on this problem of trend to the Maxwellian 
equilibrium for the fluid-kinetic models such as Navier-Stokes-Vlasov-Fokker-Planck |13j and 
Vlasov-Fokker-Planck-Euler equations [6]. 

In the next section, we introduce some notations and state the main results. 



2 Notations and Main Results 

Here are some notations which will be used throughout the paper. 

:= / fgdv, 

151^:= / {\Vvg\'^ + '^(.v)\g\'^) dv where u{v) := l + lv]"^, 
Ml ■■= / {\V^g\^ + u{v)\g\^)dxdv, 

II ||2 II ||2 II ||2 II ||2 

11511 := IIS'llL2(]g3xR3) or ||a|| := ||a||^2(K3), 
llfflloo := lbllL°°(MixR3)) 

We recall the basic properties of the linearized Fokker-Planck operator L in 1)1.50 . for instance 
see [illH]. 

{f,Lg) = {Lf,g), Ker L = Span{^}, L{vy/Ji) = -Vy/Jl. (2.1) 
We further introduce projections onto yfjl and VyJJl. 

Pi5 := {g,v^)-v^ =: u{t, x)-v^/Ji, 
P :=Po+Pi. 

Here a represents the density fluctuation and u the velocity fluctuation. Now the coercivity 
of —L can be written in terms of projections as follows (for instance see [6j): for a positive 
constant Aq > 

{g,-Lg)>Xo\{I-Po}g\l , or 
{g,-Lg)>Xo\{I-P}g\l + \u\\ 

We next define the instant energy functionals and dissipation. 

Definition 1 (Instant energy). For N > 3, there exists a constant C > such that an instant 
energy functional 8]\j{g,'SJ (p) =: ^^{t) satisfies the following: 

\a\<N 

where £N{t) will he defined in |^.i[ ). 
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Definition 2 (Dissipation). For N > 3, the dissipation rate 'D]\}(t) is defined by 
^Nit):=l (llPi^^sf + Ao||{I-PK<7||^) 

\a\<N 
\a\<N-l 

where k = min{Ao/2, 1/8}. 

Remark 3. Note that the first part of the dissipation is due to the collisions as seen in 
(j2.2p . What is distinguishable for the VPFP system is that the macroscopic part as well as 
the potential part also dissipates because of the extra damping due to Coulomb interaction. 
The exponential convergence rate is obtained through such damping of the density fluctuation, 
which can he acheived by means of the macroscopic equations. 

We are now ready to state the main results. The first theorem is about the uniform energy 
estimates for the VPFP equations. 

Theorem 4. Let N > 3. There exist an instant energy functional ^Nit) and a sufficiently 
small eo > so that if £]\f{0) < eo, then the smooth solutions {g,V(j)) to the Vlasov-Poisson- 
Fokker-Planck equations (II. 4p satisfy the following energy inequality 

j^£N{t)+VN{t)<0. (2.3) 
In particular, we have the global energy bound 

sup (Sat (t) < <fAr(0) . 

i>0 

The next theorem is about the exponential convergence to the equilibrium of the VPFP 
system. 

Theorem 5. There exist eg > and r/ > such that for small initial data ifAf(O) < eo, the 
solutions decay exponentially 

£N{t)<£N{0)e-^', (2.4) 

where rj can he chosen as 

Theorem[5l where an explicit convergence rate is given, is a direct consequence of Theorem 
Uland here we provide the proof of it. From the definition of £n ~ see (j4.ip . we first note that 

£N{t) <\y. mgf + mvct>f) . 

|a|<7V 

Next, from the Poisson equation in ()1.4p . we deduce that when \a\ = N, ||c?"V(/)p < 
IIPqC^x <?|P for |/3| = N — 1. Hence by the above definition of Djv, we see that 

Therefore, the exponential decay (j2.4p follows from the energy inequality (|2.3p by letting 
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Remark 6. The above theorems are also valid for the periodic domains without any changes 
in the proof. Our results show that the convergence rate to the equilibrium of VPFP system 
is the same exponential for both the periodic box and the whole space in the framework of 
smooth solutions. This is not the case for other kinetic models such as Boltzmann equations 
and for fluid-kinetic models [T^ . 

Remark 7. The global existence of solutions to (II. ip follows from the a priori global energy 
bound by rather standard method. In this paper, we focus on proving the uniform estimates. 

The proof of Theorem d] is based on the novel uniform- in-time energy methods developed 
for the study of Boltzmann equations over the years; for its original idea of the proof relevant 
to our model, for instance see [14j and also see [17]. It consists of the instant energy estimates 
for En and the macroscopic equations. It is interesting to see how this mechanism through the 
macroscopic equations takes full advantage of the self-consistent Poisson interaction, which 
is a key of getting the exponential convergence. We remark that the method will not give 
the exponential convergence rate for the linear Vlasov-Fokker-Planck equation in the whole 
space without the Poisson interaction. 



3 Energy estimates 

In this section, we will derive the energy estimates for the VPFP system. The first part is on 
the instant energy estimates and the second part is on recovering the full dissipation of the 
perturbation via the macroscopic equations. 

3.1 Instant energy estimates 

The goal of this subsection is to prove the following instant energy inequality: 

Proposition 8. There exists a constant C > independent of t such that the smooth solu- 
tions {g, V(/>) to the Vlasov-Poisson-Fokker-Planck equations (jl.4p satisfy the following energy 
inequality 

l^j/N{t) + VN{t) < C{£N{t)f^VN{t), (3.1) 

where 

£N{t) := Yl iP^gf + l|5?V</.f ) and 

V^it) := + Ao||{I - V}d^xg\\l) . ^'^ 

\a\<N 

Proof. First, we project ()1.4p onto {y^^} to get the conservation of the density fluctuation 

o-t + V-u = 0, -A(/) = CT (3.3) 

which will be frequently used throughout the argument. We will start with the case of a = 
in (13. ip : multiply (II. 4p by g and integrate over x M'^ to get 



1 d 
2di 



(i) (ii) 
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For (i) = f u ■ V(j)dx, we use the integration by parts and (|3.3p to get 

(i) = — / (V • u)(j)dx = / atcpdx = — A.cf)fcf)dx = I V(/)t • Vcpdx 



and thus we derive that 



For (a), we use the macroscopic variables and the microscopic part to rewrite g: 

/2 p 
—v-V(j)dxdv= i -{a^/Ji + u ■ v^/Ji + \1 — V}g)'^v -V cj)dxdv, 

then by taking the L°° of the potential, we get 

1(^01 ;^l|v</)||oo(lk||i. + Hi. + ||{i-p}5l|^). 

With (12.2p . we obtain the following: 



ld_ 
2di 



(IMP + l|V</)f) +Ao||{I-P}5||^ + IHP 



l|V.^||oo(lkf + ||nf + ||{I-P}5||2). 



(3.4) 



We next derive the higher-order estimates. Let \a\ < N. Take 5" of (II. 4p 

d^gt + V ■ V^d^g + v^-d^V,(t> + [{^g - V,g)-V,(t] - Ld^g = 0. 

The linear terms can be treated in the same way as in the case of a = 0. Hence, by multiplying 
by d^g and integrating we get 



1 d 
2di 



{Wd^gf + \\d^V<t>f) + Ao||{I - Pj^^^ll^ + 



l/3|<l«l 

= E + E =■■(!) + ill), 



dxdv 



where [7] denotes the Gauss number i.e. the greatest integer less than or equal to 7. 



Now by using the identity dx ^ g = dx ^ (y\f[i + d[ 
Cauchy- Schwartz inequality, we see that 



U ■ V 



^ + {I - Y']dx ^g and by 



\{I)\7< E l|9fV.<Al|oo(||5>f + + ||{I-PR"5||^ 



and 



m)\:< E 



+ W.-^af + W.-^uf + ||{I - V]dr^g\\l) 



+ sup|{I-PK-^g|^,) 



air + Il9>f + ||{I - PR°5||^ + V.</>f ) . 



ia„\\2 

X 
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Since > 3, Sobolev embedding yields 

E ii^f v.'Ailoo + E (ii9r''^iioo + ii5rNioo+sup|{i-p}ar^5L.);i(^7v)i 

Therefore, we obtain 

where we have used the eUiptic regularity of the Poisson equation in p.4|) 

||(/'IU»+2 < C||a||//. for s > 0. (3.5) 
This completes the proof of Proposition [HI □ 

3.2 Macroscopic equations 

In this subsection, we will prove that the macroscopic part Po5 as well as the potential part 
V0 of the solution to the VPFP system also dissipates due to the extra damping through 
the Coulomb interaction. This will be established via so-called the method of macroscopic 
equations. 

Proposition 9. There exists a constant C > independent of t such that the smooth solu- 
tions (g, V(/>) to the Vlasov-Poisson-Fokker-Planck equations ()1.4|) satisfy the following energy 
inequality 




P\<N-l 



+ C{SN{t))'^ (||afV0|p + ||afaf), 

\P\<N-l 

where 

G{t)= E /[5f^H5.,9fa + 9fa.>) + i(|5fV0p + |afap)]dx. (3.7) 

\P\<N~l 

Proof. Our starting point of the proof is the macroscopic equation for u = {g,v^/JI), which 
will be derived shortly, in addition to the macroscopic equation for a = {g, y/JI): (|3.3p . First 
we note that since {v'^^,v^ ^) = 6^^ and {v^v^ ^,v^ ^) = 0, 

= d^^a + d^^{v'v^^,{l-V]g) 

and 

{v'^, C^g - d,,g)d.,^<fy = \{v'v'^. [o^ + u^ ■ v'^ + il - P}g])d,,cP 

+ {^''VJ^ - Iv'v^VJ^, [a^ + n'^ • v''^ + {I - P}^] )d,^<P 
= {5'^^, [a^ + ■ + {I - P}g] )d,^ct> 
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Hence by further using (j2.ip . the projection of the VPFP equation (jl.4p onto {v'^^/JI} can be 
recorded as follows: 

ui + d,^a + 9,,(/> + ad.^cP + u' + 9,^ {v'v^^, {I - P}^) = (3.8) 

which is the macroscopic equation for u which is not decoupled from the microscopic part. The 
idea is to estimate a and via their elliptic coupling and the dissipation of the microscopic 
part {I — PqIs- We will prove ()3.6p first when /3 = 0. To obtain the estimate of cr = Pog 
part, we multiply (j3.8p by 8x^(7 and integrate: 



J u\dxi(Tdx + j |3s.cj|^dx + j dxi4>dxiCrdx + j adxi4>dxiCrdx + j u^'dx^adx 

+ j d., {vVVJ^, {I - P}g)dx^adx = 0. 



(3.9) 



We denote it by (/) + (//) + (///) + {IV) + (V) + {VI) = 0. We do the integration by parts 
and use p.3p to estimate (I), {III), {IV), {V) as follows: 

(/) = ^ J vJ-dxiOdx - j u^dxiCFtdx = ^ J u'dx^adx + J dxiU^atdx 

dtj y ' ' ' 

{III) = - j A^adx = j a"^ dx, 
{IV) = -\j ^^^^ 
{V) = — J dx^u^crdx = J atcrdx ~ 2^ J 
We also see that by Cauchy-Swartz inequality 



\{VI)\<^-\\Vaf + ^\\{I-P}dxg\\l 



Hence we obtain the following: 



^ I {u^dx^a + ^a'^)dx + / a^dx + ^ / \Va\'^dx 



< 



^ I \afdx + I \V ■u\^dx+^\\{I-P}dxg\\l. 



(3.10) 



In order to get the estimate of V<^ part, we now multiply p.8p by dxi4> and integrate. 
Then we get the similar expression as ()3.9p 



J u\dxi(t>dx + J dxiCrdxi(pdx + j \dxi<t)\^ dx + j a\dxi<t)\^ dx + J u'dx^cpdx 

+ J dx^ {v'v^^, {I - P}g)dxMx = 0, 

denoted by {i) + {ii) + {iii) + {iv) + {v) + {vi) = 0. The second term {ii) is the same as (///) 
which gives rise to J a^dx. The third term {iii) is a good damping term, the fourth term 
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{iv) is a nonlinear term, and we will use Cauchy-Schwartz inequality for (vi) as before. The 
fifth term (f ) forms the time integral as before: 



For the first term (i), we do the integration by parts and use (|3.3p to derive 

(0 = ^ J u'dxiCpdx - J v}dx^(t)tdx = ^ J u^d^^cpdx + J dxiU'-(j)tdx 

- !L 
~ It 



u'''dxi4>dx — / |V(/)j| dx. 



Instead of estimating the time derivatives directly, we rewrite the second term by using the 
equation again: since (j)t = —A~^at = A~^(V • u), 

j \Vc^t\^dx = j |VA-^(V-u)|2da; < j \u\'^dx, 

where we have used the boundedness of VA^^V- which is a well-known singular integral 
operator: Riesz potential [21j. Thus we get 

^ j {u'dx^<t) + ]p<i)\'')dx + j a^dx + ^ I \V<t>\^dx 

< J \a\\V^\Ux + J \u\^dx + ^\\{I-P}dxg\\l. 

By combining it with (|3.10p . we obtain 

^ J [u\dx^cl, + dx^a) + ^{\V<P\^ + a^)]dx + 2 J a^dx + ^ J {\V<P\^ + \Va\^)dx 

< j \a\{\V(t>\'^ + a^)dx + j \u\'^dx + j \V ■ u\'^dx + - Pjd^gWl. 

We next handle the higher order estimates which are needed to control nonlinear terms 
within our energy functionals. We will first derive the estimates for dx^d^a for |q| < — 1. 
Take of (|3.8p . multiply by 8x^0^(7 and integrate to get 



d"uldxid^adx + / {dx^d^ crl"^ dx + / dx^d^cpOx^d^crdx + / d'^[(7dx^4>\dxidxCrdx 

J J I (*) 

+ j dydx,dydx + J dx^d^{v'v^^,{I-F}g)dx,dydx = 0. 

We follow the same procedure as in a = case. The linear terms can be treated in the same 
way. For the nonlinear terms (*), note that 

(*) = / dydx,4>dx,dydx+ I adxX<t>dx,dydx + [ dx~^<^dxM^dx.d>dx 

= ^ j cr\d"a\^dx - j dx^crdxid"(f)dxCTdx 

^ C^iJ dx^-^aOxM^dydx- J d^-fad^^ad^dx] 



0<|/3|<|a| 
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and thus by Sobolev embedding, we get 



\{*)\ :< {£n)'^ E II^^^^I 



and in turn, 

I / {dy^d,^dy + ^\dy\')dx + J \dy\'dx + ^ J \vdy\'dx 

< J\V-dy\'dx + Ci£N)'^ E \\d^.^f + l\\{i-PW^9\\l. 



< 



l/3|<l" 

The estimates of dx^d^(j) can be obtained in a similar way: 

\dy\'dx + Ci£N)'- E (ll^^fV<Af + l|9fcTf) + i||{I-PRA"9ll^ 

m<\»\ 

By adding the above two inequahties, we finish the proof of Proposition [9l □ 

4 Proof of Theorem |4] 

We are now ready to prove Theorem U based on Proposition [8] and [H 

Proof of Theorem\^ Let k = min{Ao/2, 1/8}. By combining ()3.ip and (13. 6p . we first deduce 
that there exists a constant C > so that 

^^{8^{t) + 2KG{t))+VM{t) + 2K E (2||9faf + i||5fV<i6f + l||V9faf] 

|/3|<JV-1 ^ ^ 

< C(^jv(t))^:Pjv(t), 

where <5Ar(t) is given in (j3.2p and in (j3.7p . Since |G(t)| < £N{t), we see that |<5Ar(t) < 
^Ar(t) + 2KG(t) < ^£N{t). Now we redefine an instant energy by 

£N{t) ■.= iNit)+2KG{t). (4.1) 

Then by applying a standard continuity argument, we finally deduce (|2.3p by setting £n 
sufficiently small initially. □ 
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